Nonclassical photon pair production in a voltage-biased Josephson junction 
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We investigate electromagnetic radiation emitted by a small voltage-biased Josephson junction 
connected to a superconducting transmission line. At frequencies below the well known emission 
peak at the Josephson frequency (2eV/h) , extra radiation is triggered by quantum fluctuations in the 
electromagnetic environment. For weak tunneling couplings and typical ohmic transmission lines, the 
corresponding photon flux spectrum is symmetric around half the Josephson frequency, indicating 
that the photons are predominately created in pairs. By establishing an input-output formalism for 
the microwave field in the transmission line, we give further evidence for this nonclassical photon 
pair production, demonstrating that it violates the classical Cauchy-Schwarz inequality for two-mode 
flux cross correlations. In connection to recent experiments, we also consider a stepped transmission 
line, where resonances increase the signal-to-noise ratio. 

PACS numbers: 85.25.Cp,74.50.+r,42.50.Lc 



A voltage biased Josephson junction (JJ) produces an 
oscillating supcrcurrent, known as the ac Josephson ef- 
fect [1—3]. In a resistive environment, this current cre- 
ates electromagnetic (EM) radiation, usually in the mi- 
crowave regime. The radiation has been studied inten- 
sively in literature [3-8], not least from a metrological 
perspective, since its frequency / is given by the applied 
voltage V, in the simplest form as / = 2eV/h, where e is 
the electron charge and h is Planck's constant. For usual 
applications, involving JJs embedded in a low-ohmic en- 
vironment at moderate temperatures, a classical treat- 
ment of the interplay is adequate. Here, the microwave 
power spectrum has a peak at / broadened by thermal 
fluctuations in the bias line. At very low temperatures 
the quantum fluctuations of the resistor set a lower limit 
to the linewidth [9], through a shot noise in the charge 
transport [10]. 

In this Letter, we investigate the microwave field cre- 
ated by very small JJs, where the charge transport 
takes place through an incoherent sequence of tunneling 
Cooper-pairs [11—15]. Here, higher-order emission pro- 
cesses can also be detected, without being masked by 
classical fluctuations. This was recently addressed ex- 
perimentally [15], by simultaneous measurement of dc 
current and the power spectrum. It was found that the 
radiation consists either of a single photon or a pair of 
photons with the total energy J^. hfi = 2eV, correspond- 
ing to the energy loss of a single tunneling Cooper-pair. 
The two-photon emission is triggered by quantum fluc- 
tuations, and is of special interest, as this radiation is 
of nonclassical type. Such production of nonclassically 
correlated photons has been subject of an active study in 
various systems taking advantage of non-linearities cre- 
ated by JJs [16-23]. It has a wide interest in the fields 
of quantum communication [24] and metrology [25-27] . 

We develop a theoretical framework for perhaps the 
simplest nonclassical microwave radiation-source based 



on JJ physics, a dc voltage-biased small JJ. In partic- 
ular, we establish an input-output theory for the mi- 
crowave field in the transmission line terminated by a JJ. 
This provides a straightforward method to calculate any 
EM-field correlations in the output radiation in the case 
of weak Josephson coupling. Especially, by calculating 
photon-flux correlations in the output field, we show that 
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FIG. 1: (a) A superconducting transmission line terminated 
by a small Josephson junction. We consider an ohmic trans- 
mission line (Zo = Z\) and one with a step-like characteristic 
impedance (Z\ S> Zo), supporting modes of a A/4-resonator. 
(b) The equivalent lumped element model, (c) The emit- 
ted radiation, either single-photon (I) or two-photon emis- 
sion (II), originates from the electrostatic energy released by 
a tunneling Cooper pair, (d) The output photon-flux density 
vs drive frequency, ujd, (in units of I^Zoj%h\J^) due to classical 
Josephson radiation (blue), incoherent Cooper-pair tunneling 
(red), and thermal radiation (green, fiuj = 2Ej, Zo = 50 Q, 
T — 50 mK). The first two differ in emission at frequencies 
well below uj^. There the incoherent-Cooper-pair tunneling 
part is symmetric around the half drive Josephson frequency 
cjd/2, indicating emission due to pair production of photons. 
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the radiation is indeed of nonclassical type, as it violates 
a classical two- mode Cauchy-Schwarz inequality [28-30]. 
For optimal detection of this, we also consider a stepped 
transmission line [15], which leads to an increase in the 
signal-to-noise ratio by enhancing emission at specific fre- 
quencies. The estimated rate for creation of photon pairs 
is of the same order as in the recent two-mode squeezing 
experiments, made in the context of dynamical Casimir 
effect [18, 31, 32], showing that the detection of this non- 
classical photon production, as a side product of the ver- 
satile ac Joscphson effect, is in the reach of present tech- 
nology. 

The system we consider consists of a superconducting 
transmission line (TL) terminated by a small Josephson 
junction characterized by its capacitance Cj and criti- 
cal current I c , see Fig. la. The line is dc- voltage bi- 
ased, with a voltage smaller than the superconducting 
gap, eV < 2A. We focus on the EM radiation in the 
semi-infinite TL, which we describe by its magnetic flux 
in a discretized circuit model, Fig. lb. The TL is char- 
acterized by its capacitance Co and inductance Lq per 
unit length. We now treat the case of a homogenous TL 
and discuss the stepped impedance case at the end of the 
paper. In the continuum limit Sx — > 0, the Heisenberg 
equations of motion for the flux field operator t) has 
the form of a Klein-Gordon equation for a masslcss par- 
ticle [32], with the traveling- wave solution 
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Here Z = y/ Lq/Cq is the characteristic impedance 
and k u = u-sj CqLq the wave number. The pho- 
ton operators a! n ^ out (o;) [ai n / ut ( w )] create (annihilate) 
a photon of frequency u) moving rightwards (in) or 
leftwards (out) and satisfy the commutation relations 
Kn(out)(w),af n(out) (w')] = 6(w-u)'). The left- and right- 
moving part of the field are connected by the boundary 
condition imposed by the JJ at x — d, 
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Here <$>(x = d,t) is the magnetic flux across the junc- 
tion, $o = h/2e is the magnetic flux quantum, and 
LUd = 2eV/ h is the Josephson frequency. 

For a steady state (in the Heisenberg picture) one can 
solve for the output operators a out (ui) as a function of the 
input operators a- m {ui) . We seek a solution in powers of 
the critical current I c by multiplying the right-hand side 
of the boundary condition (2) by £ and correspondingly 
write the solution for the outgoing wave in Eq. (1) as 



flout (w) = X)^Lo £" a n(^)- The input field is independent 
of £. The zcroth-order solution (at x = d) describes the 
phase shift given by reflection at the junction capacitance 



aoM = [x(w)/x*M]a in (w), 



(3) 



where x( w ) = V^o/(l — iZ Cjuj). 

The effect of Cooper-pair tunneling appears in the 
leading-order solution 

ai (w) = ic r dt e iut sh# (t) - u d t] , (4) 
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is the the zeroth order expression for the phase fluctua- 
tion operator at the junction <f>(t) = (2n/$o)$>(d,t), giv- 
ing the dependence on the input. 
To second order in I c , we find 
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is a solution to the equation <f>\(t) = I c [4>q (t) , z] , where 
4>i(t) is the first-order result for the phase difference at 
the junction. From Eqs. (3), (4) and (6) all correlation 
functions of the output field can be calculated, to second 
order in I c . This is the main theoretical result in this 
paper. As a consistency check, we have also verified that 
these perturbative expressions for the output operators 
satisfy the canonical commutation relations. 

We can now straightforwardly calculate, e.g., the 
photon-flux density of the output radiation, defined 
as [33] 
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up to second order in I c . The result for typical 
transmission-line parameters is plotted in Fig. Id. It 
can be represented as a sum of two contributions: a 
part describing thermal radiation /th(^) and a part de- 
scribing radiation originating in incoherent Cooper-pair 
tunneling /t(w), which dominates for hui 3> k^T. We 
now introduce the tunnel impedance [13] Re[Z t (u;)] = 



|x(w)| 2 = ^o/ l + (w/w c ) 



, where the cut-off frequency 



u> c = 1 /ZqCj is given by the inverse i?C-time of the junc- 
tion and is much higher than any other relevant frequency 
in the problem. The expression for the tunnel induced 
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photon flux density now reads 
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and J(t) = ([(f>o{t) — ( / ) o(O)]0o(O)) is the phase correla- 
tion function at the junction. The function P(E) is 
the probability density for exchanging a total energy E 
with the EM environment in a single tunneling event. 
Here the term P(tvjj& — tku) describes radiation com- 
ing from forward-direction Cooper-pair tunneling and 
the term P(—huj c \ — tun) from the Cooper-pair tunnel- 
ing against the voltage bias, a process completely sup- 
pressed at temperatures fc^T < 2eV. The contribu- 
tion /t(w) agrees with Ref. 15, where the same form is 
obtained by applying the theory of incoherent Cooper- 
pair tunneling [13] to calculate the associated photon 
flux. The simultaneous electric current has then the form 
I(V) = (irhI 2 /4e)[P(2eV) - P(-2eV)}. We note that 
Eq. (9) for the photon flux is actually valid also when the 
transmission line has resonances, which then show up as 
peaks in the P(E) and Z t (uj) functions, a case treated at 
the end of this Letter. 

We now focus on the nonclassical origin of the photon 
flux for frequencies ksT < tun < fund — fcsT, where the 
simple TL environment allows for analytic expressions. 
At zero temperature, the classical result is P(E) = 5(E) 
and the radiation is created at u — m d with the total 
power I 2 Re{Z t (wd)]/2. In a quantum treatment of the 
TL, vacuum fluctuations broadens P(E) towards posi- 
tive energies and there will be photons emitted also at 
lower frequencies in < We use the long-time approx- 
imation for the phase correlations [12, 34, 35], J(t) = 
—2p [ln(u; c |i|) + 7 + z^sign(i)] , where p = Z /Rq is the 
ratio between the TL impedance and the quantum resis- 
tance Rq = h/Ae 2 and 7 is the Euler constant. Then, for 
typical transmission line parameters, p <C 1, the resulting 
finite tail at lower frequencies has the simple expression, 
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This part of the photon-flux density is symmetric around 
the half frequency cjd/2 indicating that output radia- 
tion at these frequencies occurs through pair production 
of photons symmetrically around Wd/2, Fig. lc-d. This 
symmetry is a central result of this work, and also pre- 
vails to the case of resonant transmission line. 

Photon pair production in the microwave regime can 
be achieved also through parametric effects in driven J J- 
systems, through conversion of drive photons into the 
photon pairs [16-23]. In the present case, a static volt- 
age bias is used and each photon pair is instead connected 



to a tunneling Cooper pair. We also note that the shape 
of the spectrum in Eq. (11) is inverted compared to the 
one obtained in the case of the dynamical Casimir effect 
(DCE)[31] which is oc ui(uid — uj). The difference arises 
since the Cooper-pair tunneling couples to phase fluc- 
tuations which are proportional to 1/tn, while for DCE 
the effective boundary inductance is modulated, which 
couples to current fluctuations proportional to ui. 

To verify that the photons are indeed created in 
pairs, we need to evaluate the second-order coherence 
g( 2 )(0)[28] of the out-field, which gives the probability to 
measure two photons simultaneously normalized to the 
photon flux. To second order in I 2 (and for p« 1 and 
fcgT <C hud as above) we find 

g(a)(0) = (alA ut a out a t ) ^ l + f^_\\ (12) 
« lt a out ) 2 V«Q^c/ 

which for low powers (oc I 2 ) can be made arbitrary large, 
indicating that indeed photons are emitted in pairs. 

To test whether the out-field also possess nonclassi- 
cal correlations we calculate the corresponding Cauchy- 
Schwarz (CS) inequality [36]. In a two- mode case the 
inequality reads [28, 30] 
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where the single-mode second-order coherence is defined 
as above (0) = (a\a\aiai) / ' (a\ai) 2 , the photon-flux 

cross-correlator g^(0) = {a\a a a\ab) j '{a\a a ) (a\a,b) , and 
the indices a and b refer to two different modes. We 
consider two frequency separated modes obtained by de- 
tecting the out-field in a small frequency range 8uj around 
u) a and we, [29]. On the left-hand side (LHS) we have the 
squared probability of simultaneously observing a pho- 
ton at frequency u a and wj,, which for a classical field 
is bounded by the right hand side (RHS) product of the 
probabilities of simultaneous pair detection at the indi- 
vidual frequencies. 

With the assumption ui a ,u>b 3> k^T and in the limit 
5lu — > we get to second-order in I c the equivalent CS 
inequality for the out-field 
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We notice that the nonclassicality of the field is deter- 
mined by the P(i?)-function only. Now, at zero tem- 
perature the probability to absorb energy from the EM 
environment is zero, i.e., P(E < 0) = 0. This implies 
that the RHS of Eq. (14) is zero for either ui a > u>d/2 or 
ojb > tJd/2. The LHS is finite for to a + u>b < ojd, giving a 
large regime with nonclassical correlations. For ui a = Lot 
both sides are equal, corresponding to degenerate para- 
metric down conversion, and the CS inequality is not vio- 
lated. This is the leading-order (I 2 ) result, when all cor- 
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FIG. 2: (a) Visualization of violation of the classical CS in- 
equality for flux cross-correlations. Negative values are a 
sign of nonclassicality. The photon pair production is be- 
hind the strong violation nearby the diagonal u> a + <^6 = 
Ud = 50 fj,eV/h. We plot here i>/|Max{u}|, where v — 
Sgn(d) ln(l + |d|) and d/(fieV) 2 = -P [h(w d - u a - m)f + 
P [H(tJd - 2u a )] P [h{ujd - 2u b )]. We have Z = Z t = 50 Q 
and T = 50 mK. (b) Photon-flux density due to Cooper- 
pair tunneling in the neighborhood of the optimal drive 
point cjd = ujo + w\ ~ 141 fj,eV/h, when Z\ — >• 10Zo and 
Ej — (h/2e)I c — 10 jiieV. Here, the photon pairs are mainly 
emitted at wq and wi ~ 3ujo- (c) The measured violation 
of the CS inequality d x m [m — —ft(w a )ft{wb)] in the 
diagonal uid = 0J a + of (a), (d) The measured viola- 
tion of the CS inequality corresponding to (c), in diagonal 
ujd = Ld a -\-u)b = uiQ+ui. The resonance occurs when u> a = 0Jo- 



related photons originate from a singlc-Coopcr-pair tun- 
neling event. Also at finite temperatures, the inequality 
in Eq. (14) may still be violated, as visualized in Fig. 2 a. 

We estimate the photon-flux density in the consid- 
ered ohmic transmission-line setup to be of the order 
/(wd/2) ~ 10~ 4 — 10~ 3 photons per second per band- 
width, which should be readily detectable with a state- 
of-the art experimental setup [18]. However, in the mea- 
surements of second order flux correlation functions the 
signal-to-noise ratio decreases substantially. Thus, to 
verify the nonclassicality of the field it is favorable to use 
parameters maximizing the product of the photon flux 
at the detection frequencies uj a and w&, i.e., maximiz- 
ing ft(u> a )ft(wb)- Here, we also note that the out-field 
does not possess any steady-state single mode or two- 



mode squeezing, in which case nonclassicality could have 
been shown measuring correlation functions to second or- 
der in the amplitude instead of flux[17, 19, 20, 22]. The 
phase fluctuations across the junction actually influence 
the field in much the same way as pump phase fluctu- 
ations in a parametric amplifier, making the amplitude 
correlations short-lived. 

To increase the photon flux, we consider resonantly 
shaping the P(_E)-function using a simple step struc- 
ture in the characteristic impedance [15], see Fig. 1 a. 
For Z\ 3> this provides modes at frequencies u> n = 
(l + 2n)wo. Evaluating the output field at x = 0, through 
similar theory as discussed before we obtain that the 
main results Eqs.(9) and (14) prevail, after the change 
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where C(u) = 1 + iZiCjuj and = u)\JC\L\. This 
changes the zeroth order phase fluctuation operator in 
Eq. (5), and thus P{E) through Eq. (10) and the tunnel 
impedance Re[Z t (w)] = |x(w)| 2 . 

By aligning the detection frequencies with the first two 
modes uj a — ljq and ujf, = u>i and driving at uj^ = uj a + lo\, 
the photon flux can be enhanced dramatically. In Fig. 2b 
we plot numerical results for the output photon flux den- 
sity for a setup with Z\ = IOZq, nearby the optimal 
drive for two-photon emission. The output flux is mostly 
confined into the two frequencies (of the modes), being 
even larger than the flux at the drive frequency uj^. In 
Fig. 2c-d we then plot the violation of the CS inequality 
at = uj a + U)b for the open-space configuration and the 
resonant setup, multiplied by the product of the photon 
flux at the two detection frequencies. 

In conclusion, we have derived expressions for the out- 
put field operators in a one-dimensional transmission line 
terminated by a voltage biased Josephson junction, to 
second order in the critical current of the junction. Using 
this formalism, we have confirmed the expression for the 
photon flux density derived in Ref. [15] using a different 
formalism. Furthermore, by calculating second order flux 
correlation functions we have established that the pho- 
tons below the Josephson frequency are mainly emitted 
in pairs and that the field is indeed nonclassical. Finally, 
we discuss the possibilities to enhance the photon flux 
by creating resonances in the transmission line and thus 
facilitate experimental detection. 
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